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Abstract 

W-pair production in e+e" annihilation at 500 GeV is studied using different Monte 
Carlo generators and the proposed detector for the Linear Collider in its TESLA ver- 
sion. Semileptonic W decays are used to give hints for the detector optimization and to 
determine the expected accuracy of the W mass determination and the measurement of 
the anomalous triple gauge boson couplings. With an expected luminosity of 50 fb^^ a 
W mass measurement with an accuracy of 15 MeV should be possible. The triple gauge 
boson couplings can be determined with an accuracy in the order of 10 
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1 Cross-Section of WW Pair Production and Event Rates 



The W-pair production in e+e~ annihilation is described in the Standard Model by three 
lowest order diagrams: two s channel diagrams containing the triple gauge boson couplings 
and one t channel diagram with z/g exchange, that contributes to left-handed electrons only. 
These diagrams are presented in figure [I]. 




Figure 1 : Lowest order Feynman diagrams for W-pair production in e~^e annihilation 

In the context of this study we have compared the total cross-sections obtained from 
different Monte Carlo generators in the Born approximation The Standard Model 

parameters for the cross-section calculations as used in the different Monte Carlo generators 
are summarized in table [l|. 

Table 1: Input parameters and obtained cross-section at 500 GeV for W-pair production in 
e~^e^ annihilation 



Parameters 


CompHEP 


WOPPER 


KORALW 


ERATO 


EXCALIBUR 


Z mass (GeV) 


91.187 


91.181 


91.181 


91.19 


91.19 


Z width(GeV) 


2.502 


2.501 


2.501 


2.497 


2.497 


W mass (GeV) 


80.29 


80.44 


80.44 


80.23 


80.23 


W width(GeV) 


2.094 


2.10 


2.10 


2.03 


2.03 




0.22468 


0.22995 


0.22995 


0.23103 


0.23103 


l/a 


128.0 


128.00 


128.00 


128.07 


128.07 


cross-section 


7.556 


7.530 


7.534 


7.2 


7.356 



The results of the different generators at y/s = 360, 375, 500, 800 and 1600 GeV are 
shown in figure 0. The full line represents the results from the CompHEP program. All 
cross-sections were calculated for CC03 diagrams without beamstrahlung and initial state 
radiation (ISR). No cuts were applied to the data. All cross-sections agree within a few 
percent. 
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Figure 2: Total cross-section of WW pair production 

The cross-section including initial state radiation and beamstrahlung effects is in the 
order of 8 pb at 500 GeV. Assuming 2.5 years running period at this energy, 50 fb~^ of 
luminosity will be collected. This corresponds to 400000 WW pairs. To determine triple 
gauge boson couplings from WW events in e+e" annihilation, the charge of the W has to be 
known. This is defined through the charge of the lepton. Taking into account the branching 
ratios of W~ e'v and W~ pTv, about 30000 events of each decay are expected. 
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2 Generation of W Pairs and Background Reactions 



2.1 Generation and Selection of W Pairs 

WW pairs were generated with two different generators ERATO 0] and WOPPER [Q]. The 
triple gauge boson couplings were set to the Standard Model values. Initial state radiation 
according to and beamstrahlung using CIRCE [0] for the TESLA option were taken into 
account. 

The most relevant variable for the selection efficiency of euqq and fiuqq final states is cosOi 
where 9i is the polar angle of the produced lepton. The dependence of the detector coverage 
as a function of cosOi at generator level is illustrated in fig. ||. Due to the extremely forward 
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Figure 3: Detector coverage as a function of cos9i. The sensitivity limits for the inner 
tracker and the calorimeter are indicated by arrows. 

production of the W bosons the maximum efficiency as given by the detector geometry is 
about 85 % for the calorimeter and 60% for the geometry of the inner tracker. 
34000 WW pairs of the final state fi^u qq were generated with ERATO using CC03 diagrams 
and 30000 WW pairs decaying into e~uqq as well as and 30000 WW pairs decaying into 
fx-uqq with WOPPER including CC03 and CCl 1 diagrams. 

Muons were analyzed in the inner tracker, electrons in the inner tracker , the calorimeter and 
in the luminosity detector if they hit it. 

The isolation criteria were as follows: For muons, we required the angle between the muon 
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and the closest jet to be larger than 30°. For electrons, we required less than 10 GeV of 
energy in a cone of 23°. The final cuts for the selection of WW final states with electrons 
and muons are presented in table ^ and figure ^. The selection efficiency for the electron set 
of cuts is 0.480 ± 0.003 while for the muon set it is 0.450 ± 0.003. 



Table 2: Selection cuts as used in the WW analysis 



variable 


cut num. 


value (e sel.) 


eff. 


cut num. 


value (|U sel) 


eff. 




1 


< 0.9848+lumi 


0.91 


1 


< 0.95 


0.60 


Pi 


2 


> 20 GeV 


0.90 


2 


> 20 GeV 


0.56 




3 


> 23° 


0.87 


3 


> 30° 


0.55 


Ehad 


4 


> 0.2 * ^ 


0.77 








M(lmiss) 


5 


50 - 1 10 GeV 


0.54 


5 


60 - 100 GeV 


0.46 


M(jj) 


6 


60- 100 GeV 


0.48 


6 


60- 100 GeV 


0.45 



o 
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Figure 4: Selection efficiencies as a function of the selection cuts 
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2.2 Background simulation 



We have considered two different sources of background: background from the other decay 
channels of the WW pairs and background from processes like ee — ZZ, tt, qq. All these 
processes were generated at 500 GeV using PYTHIA [jSp , including ISR and beamstrahlung. 
All background contributions taken into account are summarized in table ^. The total amount 
of these backgrounds is 2.24%. 

Table 3: Background contributions to e+e^ W^W~ Ivqq events from the different 
reactions, corresponding to the integrated luminosity of 50 fh^^ 



Reaction 


cross-section 


events 


rate 


e+e~ 


pb 


after cuts 


% 


WW — > qqqq 


3.6 






lulu 


0.8 






T~uqq 


0.6 


260 


1.8 


qq 


12.0 


3 


0.02 


tt 


0.6 


3 


0.02 


ZZ 


0.4 


50 


0.4 


all 






2.24 



2.3 Detector Model 



The detector model that we have used in our analyses was developed at DESY-IfH Zeuthen 
and is described elsewhere [§]. The following parametizations for momentum and energy 
resolution were used in the simulation program: 
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In table ^ are summarized the main parameters used in this fast simulation program. 
All generated events had to pass this program before any analysis. 

Table 4: Detector parameters used in simulation 



Detector 


Coverage 


Threshold 


Granularity 


Resolution 


\cos6\ < 


E >(GeV) 


A9" 




a 


b 


inner tracker 


0.95 


0.2 






1.5 




ECAL 


0.985 


0.1 


0.7 


0.7 


0.10 


0.01 


HCAL 


0.985 


0.3 


2 


2 


0.50 


0.04 


fi detector 


0.95 


0.2 






1.5 




LUMI detector 


0.9976 


30 


1 


1 


0.10 


0.01 



(1) 

(2) 
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3 W Resolution 



The resolution of the lepton momentum as a function of the momentum itself as well as a 
function of cos^; is shown in fig. ^. 

A cut of 20 GeV was introduced to exclude leptons from beauty and charm decays. As can 
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Figure 5: (a) lepton momentum resolution as a function of (a) the momentum and (b) cos6'/ 



be seen from this figure, the momentum resolution for leptons with a momentum larger than 
20 GeV is in the order of 2 %. This is true for the whole inner tracker (10056*; | < 0.83 ). This 
means that we have no charge mismatch for these events. 
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Events selected with these cuts were forced to have two jets. The Durham algorithm 
[ p^ was used for the jet reconstruction. Then a kinematical fit was performed requiring 
four-momentum conservation and equal masses of the two W bosons. From the selected 
events two invariant masses can be calculated: the mass of the hadronic system is the two 
jet invariant mass Mjj and M/j, is the invariant mass of the isolated lepton and the missing 
momentum. The relevant variable for the measurement of the W mass is Mjj because it is 
better determined and its errors are smaller. Fig. || shows the mass plots for the invariant 




Figure 6: W mass spectrum from invariant jet-jet mass for the electron final state 



jet-jet mass for the electron and fig. [7| for the muon final state. 

In the fitting procedure a Breit-Wigner function convoluted with a Gaussian was used as 
can be seen in the following formula: 
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Figure 7: W mass spectrum from invariant jet-jet mass for the muon final state 



ry2 

M{x,mw,Tw,(j) = / BW{y]mw,Tw)G{y - x,a)dy. (3) 

The width of the Breit-Wigner was fixed to the width of the W as used in the Monte Carlo 

generator. The a of the Gaussian represents the mass resolution of the apparatus. We ob- 
tained mass resolutions of 3.6 and 3.9 GeV for electrons and muon final states, respectively. 
The limits for the integration are chosen to be four times the experimental resolution. The 
mean value itself is shifted to lower values in the case of the hadronic mass. This is due 
to the wrong assignment of particles inside a jet, the loss of some particles outside of the 
acceptance of the detector, detector thresholds and the ISR photon. These effects can be 
corrected statistically. The error of the mean mass value represents a measure of the quality 
of the mass measurement. The minimal error of is given by the following relation: 

AMw ~ (4) 

This is true for an ideal detector and an efficiency of 100%. Replacing Vw by the mass 
resolution (4.0 GeV) and taking into account the detection efficiency of about 35 % leads 
to an expected error of the W mass of about 18 MeV using the full available statistics of 
semileptonic final states ( e^,ir^ and t^). 
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A further increase of statistics will not decrease this value significantly because of the 
dominance of the detector systematics. Better understanding of this systematics and the 
fitting procedure can improve this error slightly. 

The resolution of 6yy , cos 6i and (pi as well as cos 9q and 4>q is relevant for the determination 
of the anomalous couplings. The last four variables are defined in the W rest frame. The 
corresponding values are presented in table ^ 

Table 5: Resolutions for W production angle and its decay products angles in the W rest 
frame 



resolution of 


ji sel. 


e sel. 


cosQw 


0.0083 


0.0082 




1.20 


1.19 


COsQq 


0.026 


0.027 


0. 


2.5 


2.6 




3.0 


2.8 


cosQi 


0.050 


0.067 


<di 


2.7 


2.6 




4.0 


3.9 



"10 2250 




A©w(°) 



Figure 8: Angular resolutions for 9^ as reconstructed from the two jets 
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The angular resolutions are all in the order of few degrees. The corresponding angular 
distributions are presented in figure ^. Each distribution was fitted simultaneously with a 
gaussian to describe the W and a second order polynomial to describe the background. 
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Figure 9: Angular resolutions of the detector for lepton and jet angles 
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4 Limits on the Anomalous Couplings 



A maximum likelihood fit has been used to evaluate the sensitivity of the anomalous coupling 
measurements at 500 GeV. Monte Carlo events were produced using the ERATO generator 
and passed through the detector simulation and the analysis program. Only CC3 diagrams 
with beamstrahlung, initial state radiation and Coulomb corrections were included. The ma- 
trix element Mj for each event passing the selection criteria mentioned above was calculated 
with the ERATO generator using the 4-momenta of the lepton, neutrino and two jets after 
the constraint fit. The minimization function is given by: 



events 

-log n 



M,\ 



(5) 



Quarks are not identified so a folding has to be done taking into account the proper 
symmetrization of the matrix element squared. The measured total cross-section is denoted 
by atot- The atot was calculated using the final selection cuts for the muon decay channel. It 



will be parametrized according to ref. [|TTj] as : 

cTtot = 5*0 + 5*1 ■ a + 5*2 ■ a^. (6) 

where a denotes generically any of the anomalous couplings investigated. One parameter 
fits were performed for different anomalous couplings. All of them are expected to be zero 
in the Standard Model. The results are given in the table ^ following the notation of [|T2|]. 

Table 6: Results on accuracy of TGB couplings estimation from ERATO 



Model 


Coupling 


Value 


Ob, 




(0.0 ±2.8)x 10-3 


Ow, 




(0.0 ± 1.4)xlO-3 


Ow 




(0.0 ± 1.7) xlO-3 



The obtained values are in good agreement with the results of reference [113|]. 



5 Conclusions 



The proposed detector is well suited for W physics. 

Better equipped forward-backward regions will improve the WW results significantly. 

A conservative estimate of the reachable error of the W mass measurement is 20 MeV. 

A conservative estimate of the errors of the triple gauge boson couplings shows that 
they are in the order of 10"^. 

This will improve over the results expected from LEP2 assuming an integrated lumi- 
nosity of 500 pb^^, by a factor of 5 for the W mass and by an order of magnitude for 
the anomalous couplings. 
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